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Abstract A simple model for the dynamics of the Magellanic Stream (MS), in the 
P_l ' framework of modified gravity models is investigated. We assume that the galaxy is 

made up of baryonic matter out of context of dark matter scenario. The model we used 
here is named Modified Gravity (MOG) proposed by Moffat (2005). In order to exam- 
ine the compatibility of the overall properties of the MS under the MOG theory, the 
observations) radial velocity profile of the MS is compared with the numerical results 

likelihood analysis is performed. We also compare the results of this model with the 
Cold Dark Matter (CDM) halo model and the other alternative gravity model that 
proposed by Bekenstein (2004), so called TeVeS. We show that by selecting the appro- 
priate values for the free parameters, the MOG theory seems to be plausible to explain 
the dynamics of the MS as well as the CDM and the TeVeS models. 
PACS numbers: 05.10.-a ,05.10.Gg, 05.40.-a, 98.80. Es, 98.70.Vc 
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1 Introduction 



The asymptotic flattening of rotation curves of disk galaxies has been explained by 
invoking still undetected forms of non-baryonic dark matter 1,2. Dark matter in the 
■ form of halo, has been successfully proposed to explain the dynamics of clusters of 

galaxies, gravitational lensing and the standard model of cosmology within the frame- 
work of general relativity. Although the currently favored Cold Dark Matter (CDM) 
model has proven to be remarkably successful on large scales [3], dark matter has not 
yet been detected after several experimental efforts. Furthermore, high resolution N- 
body simulations are still in contradiction with the observations on subgalactic scales: 
the simulations predict more satellites than what is seen [H[S], and the implied spatial 
distribution of sub-halos is in contradiction with observation [6]. However, newly dis- 
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covered galaxies in SDSS and new simulations considering the environmental effect on 
the sub-halo abundance have significantly removed this discrepancy [7][8][9]. 

In order to explain the missing matter of the Universe, other alternative theories 
of gravity have been proposed. In these models, modification of the laws of gravitation 
can explain the observed asymptotically flat rotation curve of galaxies without invoking 
dark matter. One of the most famous alternative theories is the Modified Newtonian 
Dynamics (MOND) which has been introduced by Milgrom |10] , According to this 
phenomenological theory, the flat rotation curves of spiral galaxies are explained by 
modification of Newton's second law for acceleration below the characteristic scale 
ao = 1.2 x 10~ 10 ms -2 [lllU2lfT5] . Recently, TeVeSQ, a Lorentz-covariant version of 
MOND, has been presented [14]. It cannot explain the internal dynamics and merging 
of galaxy clusters such as the "Bullet Cluster" without invoking special forms of dark 
matter [TSlfTrJ] . 

Modified Gravity (MOG) is another alternative theory which has been proposed 
by Moffat (2005) 17, 18lfl9]. It is a fully relativistic theory of gravitation which is de- 
rived from a relativistic action principle involving scalar, tensor and vector fields. This 
theory explained successfully a large number of galactic rotation curves, the mass and 
thermal profiles of clusters of galaxies, the recent data of merging clusters, and the 
CMB acoustical power spectrum data [2Tll2"51l2lfll50] . 

The dynamics of the satellite galaxies and globular clusters around the Milky Way 
could be another approach to examine alternative gravity models [31133]. In order to 
discriminate, at least in principle, between CDM, MOG and MOND theories, recently 
the orbital history of Magellanic Clouds (MCs) and 3D Sun's motion in the Milky 
Way have been studied by Iorio [341135] . The consistency of the dynamics of Magellanic 
Stream ( MS, a narrow band of neutral hydrogen clouds started from the MCs and 
oriented towards the south galactic pole ) |37 , 38 in MOND theory has been investigated 
in our previous work [31]. Here we generalize that work to test MOG and compare it 
with the results obtained with the TeVeS and CDM model. 

In this paper, we use a tidal model for the formation and evolution of the MS. 
We study the dynamics of MS for different galactic models and compare the radial 
velocity with the observational data. Modified gravitational effect of galactic disk as the 
luminous part of the MW is used to study the dynamics of the MS in the context of the 
MOG and the TeVeS. The results are compared with those obtained by including dark 
matter halo under the conventional Newtonian gravity. In this work we use six different 
sets of initial conditions that have been reported in literature (Table [T}. For each set of 
initial conditions we run a large number of simulations with different galactic potentials 
and compare the outcomes to the observation. The paper is organized as follows. In 
section [2] we give a brief review on the dark halo and disk model of the MW. In section 
13] we introduce the TeVeS and the MOG models. The modeling and dynamics of the 
MS are reviewed in section U Results and discussion are given in section [S] The paper 
is concluded in section [6] 
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2 Galactic model 

We use the Kuzmin potential for galactic disk. The Newtonian potential for infinitesi- 
mally thin Kuzmin disc is given by 

$ N {R,z) = — f ~° M (1) 

where a = 4.5 kpc is the scale length and M = 1.2 x 1O 11 M0 is the mass of the disk 

In order to compare the results with the prediction of CDM theory, we add an 
axisymmetric logarithmic halo whose the gravitational potential is |49| 

* { G =~\v Q 2 \og{R c 2 +R 2 + z 2 q - 2 ), (2) 

where, R c is the core radius, Vq is the asymptotic velocity and q is the flattening 
parameter: q = 1 represents a spherical halo and q 1 gives an elliptical halo. Different 
total masses of the Milky Way in the form of halo give rise to different values of Vq. 
Therefore, there are different values of Vq in the literature: Vq = 161 kms" 1 that is 
used by Jonston et al (1999) and Law et al (2005) [52H55] , Vq = 175 kms -1 that Read 
and Moore (2005) [54] used to reproduce the tidal feature of Sagittarius, or Vq = 210 
kms -1 that Kochanek (1996)[55] obtained from an analysis of the orbital motion of 
Galactic satellites. In agreement with Helmi (2004) [55] and Ruzicka et al. (2007) 66 , 
we set Vq = 185 kms -1 and R c = 12 kpc which well reproduced the kinematics of the 
Magellanic Stream. 

We choose the logarithmic potential for several reasons. First, the relatively small 
number of input parameters of EqfJ] makes the numerical calculations faster. Secondly, 
the logarithmic potential was employed in a recent study of dwarf galaxy streams to 
investigate the MW dark matter halo [56] , and thus the application of the same formula 
allows for comparison of our results. In addition, it allows for the investigation of non 
spherical model of the Galactic halo. Finally, the flattened logarithmic halo model is 
very nearly identical to the potential of Kuzmin disk in the deep-MOND regime (see 
section [3| • 



3 MOG and TeVeS 

MOG consists of three theories of gravity called the nonsymmetric gravity theory 
(NGT), the metric-skew-tensor gravity (MSTG) theory, and the scalar-tensor- vector 
gravity (STVG). MOG has been proposed by Moffat to explain the rotation curves 
of galaxies, clusters of galaxies and cosmology without dark matter 17,18,19,20,21, 
22,23,24,25 . It was shown that a parameter- free version of STVG can be obtained 
from an action principle [26]. Good fits to astrophysical and cosmological data have 
been obtained with this more recent version of STVG. An important feature of the 
NGT, MSTG, and STVG theories is that the modified acceleration law for weak grav- 
itational fields has a Yukawa-shape force added to the Newtonian acceleration law. 
However, MOG has encountered problems in the solar system [36]. In the weak field 
approximation limit, STVG, NGT and MSTG produce similar results. From the field 
equations derived from the MOG action, one can obtain the modified Newtonian ac- 
celeration law for weak gravitational fields as 19,30 
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G(r)M(r) 

9M = ^ ( 3 ) 

G(r) = G N x { 1 + «(r) [l - e~ r / r ° (l + ^-)] } , (4) 

where G(r) is the effective gravitational coupling constant, Gjv is the Newtonian gravi- 
tational constant, M is the baryonic mass, and a(r) = \J MqM~ x . Because of the large 
galactocentric distance of LMC, we use the point mass approximation for the total stel- 
lar mass (disk and bulge). The masses of the disk and the bulge are M^g^ = IO^Mq, 
M bu i ge = 3.4 x 1O 1O M , yielding a total baryoinc mass of M = 1.34 x 10 n M Q [53] . 
Moreover, Mo(r) and ro(r) are two scaling parameters that vary with distance and 
determine the coupling strength of the vector field to baryonic matter and the range 
of the force, respectively. These parameters are determined by the equations of motion 
for effective scalar fields derived from an action principle [19] . In order to calculate the 
MOG dynamics, we have to phenomenologically obtain Mq and r$: this determines 
G(r). It is postulated that Mq and tq give the magnitude of the constant acceleration 
as 

GM 

90 = — —■ (5) 

r o 

We assume that, for galaxies and galaxy clusters this acceleration is determined by 
go = cHq, where Hq = WOh kms^ 1 Mpc^ 1 is the current measured Hubble constant 
and h = (0.71 ± 0.07). This gives g = 6.90 x 10~ 10 ms -2 . 

Analyzing the galaxy rotation curves, a satisfactory fit to LSB and HSB galaxy 
data is obtained with the parameters Mq = 9.6 x 10 1:l Mq and ro = 13.9 kpc, whereas 
in the dwarf galaxies smaller than 12 kpc, the best fit value of parameters obtained as 
Mq = 2.4 x 10 11 Mq and ro = 9.7 kpc, and for the satellite galaxies, the parameters 
are M = 4.6 x 10 13 M Q and r = 111.3 kpc [2S]. An empirical fitting of M versus 
ro for the wide range of spherically symmetric systems, from the solar system scale to 
clusters of galaxies has been obtained and depicted in Fig. 2 of [30]. The modifications 
to gravity in Eq. ([3]) would be canceled by decreasing Mq and increasing ro (i.g. Mq — > 
and ro — ► oo). These parameters are scale dependent; thus, they are not to be taken 
as universal constants. MOG is not arised from a classical modification, but from the 
equations of motion of a relativistic modification to general relativity. 

According to TeVeS model, the physical metric near a quasi-static galaxy is given 
by the same metric as in general relativity with just a little change: the Newtonian 
potential for known matter, (j>N replaced by the total # which comes from two parts, 

<P = $ N + (f>s, (6) 

where, <j> s is a potential due to a scalar field. The added scalar field plays the role of 
the dark matter gravitational potential and the corresponding modified acceleration is 
9 — 9N + 9s- The scalar field potential for the Kuzmin disk is as following [51] 

4> s ^ iMGa / /2 HR 2 + (\z\ + af). (7) 

From Eqs. |l]), ([6]) and (0, we calculate the accelerations of a test object orbiting the 
galaxy. 

Because the Poisson's equation is non-linear in MOND, the strong equivalence 
principle is violated |TT] , and consequently the internal properties and the morphology 
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of a stellar system are affected both by the internal and external fields. This effect 
is specific for MOND; it significantly affects non-isolated systems and, in principle, it 
should be taken into account. 

fn the presence of an external field, the total acceleration, which is the sum of the 
internal a, and external a e accelerations, satisfies the modified Poisson equation [TT] 

V.[ M ( |ai + ae| )(a» + a e )] ~ 4ivG P , (8) 

where a e approximately is constant, a^ = V0 is the non-external part of the potential, 
and p is the density of the star cluster. For spherical system one can approximately 
write equation [8] as aj/x(|aj + a e |/ao) = a^r- The EFE is indeed a phenomenological 
requirement of MOND and it was postulated by Milgrom (f983) to explain the dy- 
namical properties of nearby open clusters in MW. Equation[S]is only an approximate 
way to take into account the external field effect, in order to avoid from solving the 
modified Poisson equation with an external source term p ex t on the right-hand side. 
The EFE allows high velocity stars to escape from the potential of the Milky Way [581 
I59| . and implies that rotation curves of spiral galaxies should fall where the internal 
acceleration becomes equal to the external acceleration [60 , 61 . 

In the three body problem of LMC-SMC-MW interaction, the external field of MW 
plays an important role in the internal interaction of LMC and SMC as MONDian EFE 
( for more details see the paper by Iorio [33)- Since the used model for magellanic 
system in this paper is the single cloud model (i.e. we neglect the presence of SMC), 
there is no mutual interaction between Clouds, and therefore we do not consider the 
EFE in the orbital motion of LMC. In the next section using the RK4 technique we 
solve numerically the equations of motion. 



4 Modeling and Dynamics of Magellanic System 

The Magellanic Stream is due to past LMC-SMC-MW interaction; it extends along a 
curved path as a narrow band of neutral hydrogen clouds, originating from the MCs 
and oriented towards the south galactic pole [371138] . The radial velocity and the column 
density of this structure has been measured by many groups [391140] . The observed HI 
radial velocity profile is measured along the Magellanic Stream in the Galactic Standard 
of Rest (GSR) frame (see caption of Table 1). Large-area 21-cm0 radio surveys have 
produced most of the information now available about the detailed structure of the MS. 
Since its discovery as a long stream of HI gas trailing the MCs, a number of models 
have attempted to explain the dynamics and origin of the MS. Observations show that 
the radial velocity of the MS with respect to the Galactic center changes from kms -1 
at the beginning to —200 kms -1 at the end of the tail. 

There are two main explanation for the origin of the MS: the tidal-stripping 42,43 
44 45 62 ,63 and the ram-pressure model 46,47,48,64,65,66 . According to the tidal 
hypothesis, the interaction of the MCs with the MW turns out the materials to form the 
tidal tails emanating from opposite side of the Clouds. In the tidal scenario, although 
the observed radial velocity profile of the stream has been modeled remarkably well, 

2 The electron and proton in the hydrogen atom can have their spins parallel or antiparallcl. 
A transition between these two states is called a spin-flip transition and leads to the emission of 
a photon whose wavelength is 21 cm. This is in the radio part of the electromagnetic spectrum. 
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IC 


3D v(x,y, z)(km/s) 


r(x, y, z)(kpc) 


GN96 
HR94 
K2 
Kl Mean 
M05 
vdM02 


(-5,-226,194) 
(-10.06,-287.09,229.73) 
(-91,-250,220) 
(-86 ± 12, -268 ± 11, 252 ± 16) 

(-4.3,-182.45,169.8) 
(-56 ± 36, -219 ± 23, 186 ± 35) 


(-1.0,-40.8,-26.8) 
(-0.8,-40.8,-27.9) 
(-0.8,-41.5,-26.9) 
(-0.8,-41.5,-26.9) 

(0,-43.9,-25.1) 
(-0.8,-41.5,-26.9) 



Table 1 The LMC initial conditions (IC) in Galactocentric coordinates which is a right- 
handed Cartesian system with the origin at the galactic center, the galactic plane at Z = 
and the sun at the position Rq= (-7.9,0,0). References: Gardiner & Noguchi (1996, GN96), 
Heller & Rohlfs (1994, HR94), Kallivayalil et al. (2006a,b, Kl and K2), Mastropietro et al. 
(2005, M05), and van der Marel et al (2002, vdM02). 



the smooth HI column density distribution does not agree with the observations and 
the expected stars in the stream have not been observed yet. 

The other hypothesis invokes the idea of ram pressure stripping of gas from the 
MCs by an extended halo of diffused gas around the Galaxy. The drag force on the 
gas between the MCs causes weakly bound material to escape and form a trailing 
gaseous stream. This material after escaping from the MCs falls towards the Galaxy. 
On the other hand, also the ram pressure models have their own problems. The clumpy 
structure of the Stream can hardly be reproduced by the process of continuous ram 
pressure stripping and can not reproduce the observed slope of the radial velocity 
profile along the MS, especially the high negative velocity tip of the Stream. Naturally, 
if the gas in the galactic halo has a clumpy distribution, which seems likely, this would 
lead to a clumpy MS in the ram pressure scenario. 

In this section we calculate the dynamics of MCs in modified gravity theories, using 
their its present positions and velocities as the initial conditions for the equations of 
motion. The overall external force on the MCs is the sum of the gravitational pull by 
the Galactic halo and disk, hydrodynamical drag force from the extended gaseous halo 
and dynamical friction force from the Galactic halo. In MOG and TeVeS models, since 
we have ignored the Galactic halo, there is no dynamical friction and the main factor 
in the gas stripping is the tidal force exerted by the Galactic disk. It should be noted 
that, in the absence of extended dark halo, dynamical friction produces if an object 
moves through the visible mass distribution of the host galaxy. There is no dynamical 
friction for an object that moves outside the mass distribution. In the case of MCs, 
since the Clouds move very far from the center of the MW (i.e. outside the visible mass 
distribution), there is no dynamical friction in the equations of motion. The general 
equation of motion of the MCs is 

d 2 r c d , , fc + Fdrag 

— = —te G (r c )]+ ^ (9) 

where, r c is the distance of the MCs from the Galactic center, m c is the mass of the 
Cloud, 4>q is the gravitational potential of the Galaxy, f c is the dynamical friction 
force, and Fdrag is the hydrodynamical drag force. We adopt the standard form of 
dynamical friction as follows [S7] 

f c = 0.428 In A^$-^, (10) 
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where In A ~ 3 is the coulomb logarithm and v c is the relative speed with respect to 
the gaseous halo. Here we model the MCs with a dense sphere moving through the 
gaseous halo; thus the drag force on this sphere is given by 

Fdrag = 'dpgV 2 c Tl D 2 , (11) 

where is the coefficient of drag force, p g is the density of halo and D is the size of 
the MCs. The drag force plays no significant role in the orbital motion of LMC [32| . 
Therefore, we will neglect the drag force for the rest of the paper. 

One of the problems with modeling the MCs-MW interaction is an extended high 
dimensional parameter space. In order to reduce the parameter space, one can neglect 
the influence of the SMC in the LMC-SMC-MW interaction (single cloud model). Lin & 
Lynden-Bell [67] showed that such configuration could explain the existence of trailing 
tidal stream. Furthermore, Sofue 48 considered the continuous ram pressure stripping 
in simulation of the Magellanic system, ignoring the presence of the SMC. Recently, the 
LMC-MW interaction without including the action of SMC has been studied, and the 
Stream's properties have been successfully reproduced [65] . The assumption of ignoring 
the SMC is also motivated by the fact that the SMC's impact on the orbital history of 
the LMC is minimal, because of AIlmc 3> M^mc EH- Thus, the global dynamics of 
the LMC alone is sufficient for discrimination of gravity model. Since the aim of this 
study is examining the overall properties of the MS in the context of MOG and TeVeS, 
we use the single cloud model and ignore the presence of SMC. 

In the modified gravity models, since we have ignored the Galactic halo, there is 
no significant dynamical friction on the MS, and the main factor in the stripping of 
the hydrogen from the MCs is the tidal force exerted by the Galactic disk. 

Another simplifying assumption is the elongation of the MS. Although the forma- 
tion of the MS is still a subject of debate, both the ram pressure and the tidal stripping 
models predict that the orbits of the Clouds should trace the MS at least for less than 
a Gyr in the past. In other words, the MS, as seen on the sky, is following the MCs and 
it is reasonable to assume that the MS is a narrow and long column of gas which is 
moving along the same orbit as the MCs ,44,45,65,69]. This assumption is motivated 
by the work by Johnston et al. [70] , who suggested that the tidal stream acts as "fossil 
record" of the recent orbital history of their progenitors and could thereby provide a 
probe of the galactic potential. In this way, the mean velocity of gas in the MS will be 
the orbital velocity of the MCs. In order to explain this assumption one can consider 
the dynamics of nearby satellites around the MW under the external tidal field 71 , 72 
73 , 74 , 75 . There is a connection among the tidal tails and satellites's orbits. After the 
formation and evolution of tidal tails around the satellite, the tails are aligned with the 
orbital path only when the satellite is near the perigalaction of the orbits |73II74| . The 
degree of elongation of the tails along the satellite orbital path strongly depends on 
the ellipticity of orbit. In the circular orbit case, the tails are clear tracer of the cluster 
path, but for most eccentric orbits, the tails are strictly elongated along the orbital 
path only when the satellite is near the perigalaction, whereas at the apogalaction they 
tend to deviate from the satellite path [7S]. Differently stated, in the case of eccentric 
orbits, the angle between tail and orbit decreases, when approaching the pericenter 
and reaches a minimum at the pericenter, and then increases moving away from it. 
Since the MCs are at the perigalaction, we assume that the stream is approximately 
elongated along the orbital path, and thus the mean velocity of gas would be the orbital 
velocity of the MCs. It must be noted that if the gas in the galactic halo itself moves, 
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Fig. 1 The radial velocity versus the angular distance along the MS plotted for differen galactic 
models and compared with the observational data. The results of the Newtonian gravity model 
without dark matter, the CDM halo for different values of q, and the best fit TeVeS model 
are indicated. The Moffat's modified gravity predictions in terms of the standard values of Mo 
and ro (M-MOG) and with the best fit values (M-MOG-bf) are also presented. Each figure 
corresponds to one of the initial conditions listed in Table [T] 



this effect can put the MS into almost any direction, making it impossible to use its 
dynamics to test CDM versus different gravity theories. 

In the next section, we will use tidal scenario for the dynamics of MS obtaining 
the trajectories of the MCs in various gravitational models and comparing the radial 
velocity profile of MS with the observed data. 
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MOG 




r ( kpc ) 

Fig. 2 Marginalized likelihood functions of the free parameters of MOG theory (ro,Mo) for 
various initial conditions. The location of picks correspond to the best fit value of parameters. 
The intersection of the curves with the horizontal solid and dashed line give the bound at la 
and 2cr confidence levels, respectively. 



5 Results and discussion 



The results of maximum likelihood analysis of the radial velocity profile of the MS for 
various gravity models, are presented in this section. One of the challenging points in 
calculating the orbital motion of the MCs is the initial conditions. The trajectory of 
the MCs is highly sensitive to the initial conditions, and slight variation may lead to 
completely different scenarios [34]. Furthermore, the parameter space of the Magellanic 
system evolution is very large [41,42,43,44,45 . Due to the implied numerical burden, 
doing the fully self-consistent simulation for each set of initial conditions and differ- 
ent combinations of parameters is impossible. Recently, an extended analysis of the 
parameter space for the interaction of the Magellanic system with the MW has been 
done [66] . They have used the genetic search algorithm and an approximate restricted 
N-body simulation method. Also in another orbital analysis, different sets of initial 
conditions have been used to determine the dynamics of the LMC |69j . 



10 



Haghi and Rahvar 




Fig. 3 Marginalized likelihood functions of the critical acceleration in TeVeS, for various set 
of initial conditions. The location of the picks correspond to the best fit values of parameters. 
The intersection of the curves with the horizontal dashed and solid line give the bounds at lc 
and 2<r confidence levels, respectively. The currently accepted value for ao = 1.2 on galactic 
scales are indicated by vertical dotted line. The predicted value of ao in the case of HR94 is 
in good agreement with the conventional value ao- Except GN96 and M05, the other initial 
condition predicts the acceptable value for ao at 2cr confidence levels. 



We numerically integrate the equations of motion in 3D using the forth-order 
Runge-Kutta technique with a timestep of 0.5 Myr and total time of 5 Gyr. In MOG 
and TeVeS models, the radial velocity profile of the MS is extracted numerically for 
various initial conditions of the MCs listed in TableQ] In the next step we use x 2 fitting 
to constrain the parameters of the model. Figure Q] compares the observed data with 
the theoretical predictions for various initial conditions. The data points represent the 
angular variation of observed radial velocity with respect to an observer located at the 
center of Galaxy |40) . The zero angular distance corresponds to the position of MCs 
center of mass. The models reproduce radial velocity of the stream as an almost linear 
function of angular distance with the high velocity tip reaching 200 kms -1 . 

To test the consistency of MOG and TeVeS with observations, we compare their 
theoretical predictions to the data directly obtained from observations and find the 
best-fitting parameters. Since the agreement with the data cannot be perfect, we give 
confidence intervals for the free parameters of the model using likelihood analysis. We 
compute the quality of the fitting through the least-squares fitting quantity \ 2 defined 
by 

2, 1 (VtheoryiPa} ~ Vobs) 2 

x{p^ = n1^ -2 • ( 12 ) 

i=l 1 

where cr, is the observational uncertainty in the radial velocity, N is the number 
of degrees of freedom and p a is the model parameters. In the case of MOG, p a = 



3 N = n — n Pa , where n is the number of observational data points and n Pa is the number 
of free parameters. 
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X 2 (9 = 0.75) 


22.9 


2.39 


2.2 


8.11 


87 


38 


9 


x 2 (q = i) 


11.2 


5.58 


2.36 


14.2 


52 


9.8 


10 


X 2 (l = 1-25) 


4.9 


10.1 


4.45 


20.7 


32 


4.24 


11 


g 


1.7 


0.4 


0.7 


0.2 


2.7 


1.43 


12 




1.61 


1.67 


2.0 


3.1 


2.83 
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Table 2 The best fit parameters of gravitational model and corresponding minimum values 
of x 2 f° r different initial conditions (IC). The second row shows the \ 2 of the Newtonian 
gravity. The \ 2 an d corresponding best fit parameters in MOG represented in rows of 3-5 
as well as TeVeS in rows of 6 and 7. The error bars are obtained by marginalized likelihood 
analysis (Figures 2 and 3). The results of logarithmic CDM halo model for different values of 
q are shown in rows of 8-10. The last two rows indicate the best fit value of q and \ 2 f° r the 
logarithmic halo. 



{ro,Mo} and in the case of TeVeS, p a = clq. Using Eq. (|12|) we find the best fitted 
values of the model parameters that minimize X 2 {??«}• 

To constrain the parameters of the model, we carry out statistical analysis using 
the marginalized likelihood function defined by 

<{pa} = £e-x 2 ^}/ 2 , (13) 

where £ is a normalization factor. The marginalized likelihood function of the model 
parameters are plotted in Figures 2 and 3. According to Eq. (|13|l . the maximum values 
of £ correspond to the minimum values of \ 2 ■ Therefore, the location of the peaks 
correspond to the best fitted values of the parameters. The main benefit of using the 
£ function is in finding the la (68.3 %) and 2a (95.4 %) confidence intervals of best- 
fitting values of the parameters. The best-fitting values for the parameters of the model 
at la and 2a confidence intervals for various initial conditions are shown in Figures 
2 and 3. The best-fitting values for the parameters of the model at the la confidence 
interval with the corresponding x 2 are presented in Table 2. 

The second row of Table [5] shows the x 2 of Newtonian gravity considering the only 
baryonic matter for Galaxy. As we expect, the values of x 2 are very large which means 
that we need the dark matter or modification of gravity. The most interesting result 
is that the initial condition M05 is almost compatible with the observation due to the 
low value of x 2 > which predicts very small amount of dark matter for the galaxy. 

In MOG model we apply the likelihood analysis to find the best values of the mass 
scale Mq and range parameter ro- The minimum values of x 2 and corresponding best fit 
parameters are represented in the rows of 3-5. The error bars obtained by marginalized 
likelihood analysis (Fig. [5]). According to the values of AIq and ro in section[3] which were 
obtained from fits to satellite galaxies and rotation curves of galaxies, the reasonable 
ranges for Magellanic system would be in the range of [90 - 5000] x 1O 1O M for M 
and [14 — 111] kpc for ro- However, the obtained best fit values of ro are almost in 
the acceptable range for all sets of initial conditions, but the Mq values are out of 
range in the case of GN96 and vdM02 (see Table [2]). Furthermore, the values of ro and 
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Mq in Kl-Mean and M05 do not follow Eq. |[5j) which postulated in MOG theory even 
considering the error bars. The case of HR94 is in excellent agreement with observation, 
because of the reasonable best fit parameters and the smaller value of x 2 - Therefore, 
we can conclude that, the MOG theory could be compatible with the observational 
feature of MS by choosing the appropriate values for free parameters and for special 
sets of initial conditions of the MCs. 

In the TeVeS model, with the standard value of the acceleration scale, ao — 1-2 x 
10 -10 ms -2 , except HR94 which reached a low x 2 value, the other initial conditions 
are not compatible with observations. However, the value of ao has been fixed from 
rotation curve analysis by Begemann et al. (1991)[T2"]. it is worth to obtain the value 
of ao from some other independent method. Moreover, among the MOND community 
there is no common idea about the value of ao at different scales. At different scales 
like clusters of galaxies or small dwarf galaxies, it seems the standard value of ao does 
not give acceptable results. At large scales, a lower value of ao, and at subgaactic scales 
a larger value of ao work better 79, 80, 81 , 82 ). In this paper, since we are studying the 
dynamics of MCs, i.e. local group scale, it is probable to find another best-fit value for 
ao- For this reason, in order to see the effect of different choice of ao, we allow ao to 
changes and find a best-fit ao for each initial condition which gives the better fit with 
observation. 

Figure depicts the best fit values of ao for different initial conditions. Except 
models GN96 and M05 which prefer a small value of ao, for the other initial conditions 
(Table |S), the best fit value of ao are almost in agreement with standard value. Thus, 
in the TeVeS model, GN96 and M05 are ruled out due to their strange perdition for 
ao. 

In order to compare the results with the CDM models, we apply the same anal- 
ysis for the logarithmic halo model. We adopted different values for halo flattening 
parameter to deal with prolate, oblate and spherical halos. According to Table [2j the 
minimum values of x 2 belong to HR94 and K2. They prefer the oblate and spherical 
halo models (q — 0.75 — 1). According to the last row, HR94, K2, and Kl-Mean prefer 
oblate halos whereas, GN96, M05 and vdM02 prefer the prolate ones. In addition, the 
minimum values of x 2 , in the CDM model for GN96 and HR94, are comparable with 
MOG and TeVeS, which means that the alternative gravity models can successfully 
explained the observational velocity profile as well as the CDM halo models. 

According to the minimum values of x 2 for different gravity models in Table 2, for 
the majority of initial conditions, the MOG model fit the observations better. However, 
the difference is not enough to discriminate between different gravity models. 



6 Conclusion 

In our paper we test alternative gravity models by comparing of the radial velocity 
profile of the Magellanic Stream (MS). We numerically integrated the orbits and dy- 
namics of the Magellanic Clouds (MCs) within different gravitational models. The MS 
is produced by the interaction of the MCs with the Galaxy, and is considered to be 
a tracer of the MCs. The drag force plays no significant role in the orbital motion of 
the Clouds. In the absence of the dark halo of the Milky Way, since the Clouds move 
outside the visible mass distribution, there is no dynamical friction in the equations 
of motion. We used the single cloud model and ignore the presence of the SMC. The 
gravitational tidal effect of Galaxy is assumed as the origin of the MS and it follow the 
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same orbit with the same dynamics as the LMC. The radial velocity of this structure is 
compared with the observation, allowing us to put constraints on the free parameters 
of models. 

A preliminary inspection of the fits to the Magellanic Stream showed that for some 
individual initial conditions, the MOG theory, choosing the appropriate values of free 
parameters, tq and Mq, can be consistent with the observational velocity feature of MS, 
as well as TeVeS, and CDM hypothesis. However, due to small difference in minimum 
values of x 2 > the discrimination between different gravity models is impractical. The 
fits to the Magellanic Stream based on the parameters Mq and tq of the older study of 
STVG, reveal that the new parameter-free version of STVG [26] appears to agree well 
with the results presented in this paper. 

We should point out that while the dynamics of the MS is in a good agreement 
with observations, a problem with tidal scenario could be the lack of corresponding 
stellar tidal debris in the MS. In TeVeS and MOG, since the tidal radius is larger than 
of the Newtonian case, we expect a Galactic structure with a gas distribution that is 
more extended than the stellar component form of the MS. N-body simulations of tidal 
stripping of the MS from the MCs in MOG without dark matter will give a better view 
of this model and enable us to compare the density distribution of the gas in the MS 
with observations. As an additional result it seems that the initial condition HR94, is 
in a good agreement with observation in all gravitational models. 

Acknowledgements We thank Alireza Moradi for useful comments on the English of paper. 
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